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Solid solutions of KTaOa and KNbOs (KTN) exhibit 
unusual dielectric properties at ferroelectric phase transi- 
tions, which are usually associated with a Nb instability 
Q that is often explained by the large Nb dynamical 
charge in oxides. Using non-self-consistent VASP com- 
putations 0, we showed earlier that the Nb dynamical 
charge in layerd KTN strongly depends on the local en- 
vironment of Nb. In the present study, we have per- 
formed self-consistent first-principles computations for 
fixed composition KTay/gNbi/gOs in different configu- 
rations shown in Fig. 1. We computed the equilibrium 
ionic positions. Born charges, force constants, vibration 
frequencies and modes. 

In agreement with our previous results, we found that 
the Nb dynamical charge is sensitive to the nearest neigh- 
bors (Table 1). The Nb charge in directions where there 
are only Nb ions is about 9.9e. In the direction where 
there are Ta ions, the Nb dynamical charge is about 8.3e. 
In contrast, the Ta charge increases in directions where 
there are Nb ions. The comparison of the energies (Ta- 
ble 2) shows that there is tendency to Nb ordering in 
planes if one does not take into account external stress 
necessary for obtaining experimental lattice constant in 
KTaOa. When we applied this stress, we obtained that 
the energies of the considered structures are close to each 
other. Note that this situation is very different from that 
known for PST, PMN etc.: in those cases, 1:1 or 1:2 or- 
der appears as a result of electrostatic interactions due 
to different charges of Nb and Mg or Sc and Ta 
KTN, the valence of Nb and Ta is the same. 



TABLE I: Ta and Nb dynamical charges in KTN 181. The 
numbers correspond to the Ta ions from the left to the right 
(Fig. 1). 
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FIG. 1: The geometry of the structures considered. 



TABLE II: Total energies (with respect to the highest value 
obtained) for the computed supercells [eV]. "equilibr." is the 
self-consistent energy obtained at zero stress. "Vegard" is the 
energy obtained at the volume which is the average between 
KTaOa and KNbOs. "KT-stress" is the energy obtained at 
the hydrostatic stress giving experimental lattice constant for 
KTaOa. 
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We obtained that the average Nb-0 distance is larger 
than the average Ta-0 distance. It implies that the Nb 
ions cause tensile distortion in KTN. Nevertheless, the 
lowest soft TO modes are mostly due to Nb (Table 3). 

We also found that the lattice instability in the KTN 
181 configuration appears in the directions of the Nb 
planes. That is consistent with other first-principles 
computations The instabilities can be explained by 



the long-range dipole-dipole interaction supplemented by 
the covalent (pseudo vibronic) effect: the displacement 
of one of the Nb ions along the Nb-O-Nb bond triggers 
appearing an additional force on nearest Nb ions in the 
same chain due to the change of electronic structure that 
is produced by the shift of the Nb ion 3|: ecor{Ri, R2) = 



E [W{Rl)r,r'W{R2)r'A F (l - F') I 
t,t' ^ 



where 



2 



TABLE III: Vibration frequencies obtained for TOl soft 
modes [cm~^] at the external stress corresponding to an ex- 
perimental lattice constant in KTaOa. 
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W(TC) = r idH / dR) is the linear vibronic operator, r 
being the ionic displacement and H being the electron 
Hamiltonian, e'^, /"^ are electron energy and occupation 
number, respectively, r denotes the electronic state. 

By using a tight-binding approach 0| we have com- 
puted the correlation energy for the ionic displacements 
in KNbOs, in {z^ / R^^ q) eV units, where z being the 
ionic charge and R — RMb-o is the Nb-0 bond length. 
When the nearest Nb ions are shifted, the correlation en- 
ergy is -10.3 and -0.9 for the displacements along and 
perpendicular to the Nb-O-Nb line, respectively. If Nb 
and nearest O are shifted then the correlation energy is 
25.9 and 8.6 for the displacemnts along and perpendicu- 
lar to the Nb-0 bond, respectively. Hence, the vibronic 
effect facilitates the displacement of the Nb ions in the 
chains in the same direction while the Nb and O ions 
are easily shifted in opposite directions. This result is in 
qualitative agreement with straightforward computations 
of the interatomic force constants (see and references 
therein), which showed that the M-0 force constant along 
the . . . -M-0-M-0-. . . chains proves to be practically neg- 
ligible while the M-M force constant appears to be sur- 
prisingly strong. 

Hence, the local polarization in KTN appears in the 
Nb-reached regions. The polarization direction follows 
the field and the free energy can be written as: F = 
a{P-Po)^+A{P-Pof + ...-{E + Eo)P where Eq is 
random and P is polarization magnitude (P > Pq > 
where Pq is the lowest possible, remnant, magnitude 
of polarization). An unusual field dependence of suscepti- 
bility follows from this: 1 /%(£') = 2a + 3AE/a + E > 
Eq. At smaller fields, susceptibility is constant or de- 
creases because of averaging over Eq . 

The strongly angle dependent correlation of the dis- 
placements can explain the diffuse X-ray scattering in 
KNbOs 1^ without the assumption that the Nb ions have 
deep off-center wells. At high temperatures (where the 
soft mode is stabilized), acoustic excitations of the (Nb 
or O) chains can be intensive due to nearly flat transverse 
acoustic phonon energy surface in (001) planes 10] . The 
peaks in the XANES data can be understood if the 
large chain-like fluctuations are considered together with 
the strong dependence of the X-Ray absorption matrix 
element on the ionic displacements |l2|. 



The dynamical nature of the cooperative displace- 
ments |l3j is supported by the fact that, in the high tem- 
perature cubic phase, there is little dielectric dispersion 
[lj|. In the ferroelectric phases, the situation is more 
complicated. For instance, in the tetragonal phase, the 
average polarization is directed along the [001] direction. 
However, the direction of a local Nb displacement is [111] 
due to anharmonicity. The correlation of the displace- 
ments can result in the appearance of finite size chains 
directed in the [001] directions, while the metal ions are 
displaced in [111] directions (common within a flnite-size 
chain but random in different chains [13 ). One can in- 
troduce population numbers for such chains correspond- 
ing to different directions of the displacements. Note that 
these (quasistatic) chains can exist only because of i) an- 
harmonicity of the potential relief, where the direction of 
the (local) displacement does not coincide with the direc- 
tion of the polarization, ii) strong correlation of the dis- 
placements, iii) existence of a uniform macroscopic field. 
In the cubic phase, the last condition fails and the disper- 
sion connected with the chains disappears. In PbTiOs, 
the low-temperature phase is tetragonal: the ions are 
shifted in the [001] direction and the first condition fails; 
the dispersion is also absent in accordance with exper- 
iment [ij]. In the ferroelectric phases of BaTiOa and 
KNbOa, all three conditions are satisfied and the disper- 
sion exists [l^. 
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